library(leaflet)
library (raster)
library(rgdal)
library(RColorBrewer)
library(devtools)
#devtools::install_github("rstudio/leaflet")
setwd("C:/Users/crish/OneDrive/Documents/Columbia SIPA/2017 Spring/Data Visualization/Final Project/")
library(readr)
iycf <- readstata13::read.dta13("data/MAP_2015_General_Info_IYCF.dta")
sam <- readstata13::read.dta13("data/2015_SAM_TOTAL.dta")
mnp <- readstata13::read.dta13("data/2015_MNP.dta")
Infant and Young Child Feeding (IYCF)
library(raster)
world_spdf <- shapefile("data/world_map/TM_WORLD_BORDERS_SIMPL-0.3.shp")
# Merge the QoG data to Shape Files
library(dplyr)
IYCF2 <- world_spdf@data %>% left_join(iycf, by = c(FIPS = "fips"))
world_spdf@data <- IYCF2
Make map…
leaflet(world_spdf) %>%
setView(lat=10, lng=0 , zoom=2) %>%
addPolygons(
stroke = FALSE, fillOpacity = 0.5, smoothFactor = 0.5,
color = ~colorNumeric("RdYlGn", Health_facilities)(Health_facilities),
popup = paste("Country:",world_spdf$NAME,"<br/>",
"Health facilities:",round(world_spdf$Health_facilities,1),"","<br/>",
"Community health workers:",round(world_spdf$Community_health_workers,0),"")) %>%
addLegend("bottomright",
pal = colorNumeric("RdYlGn", world_spdf$Health_facilities), values = ~Health_facilities,
title = "Number of health care facilities", opacity = 0.5)
Micronutrient powders
Children given Micronutrient Powder, per region, 2015
childrenMNP <- c(203552,2285242,570648,4263234,446194,1785867,612016)
Region <- c("CEE/CIS", "East Asia and Pacific", "Eastern and Southern Africa",
"Latin America and the Caribbean", "Middle East and North Africa",
"South Asia", "West and Central Africa")
MNP <- data.frame(childrenMNP,Region)
names(MNP) <- c("childrenMNP","Region")
library(ggplot2)
library(ggthemes)
library(scales)
ggplot(data=MNP, aes(x=Region, y=childrenMNP, fill=Region)) + geom_bar(stat = "identity") +
scale_fill_brewer(palette="YlGnBu") +
ggtitle("Graph 1. Number of children given Micronutrient Powder, per region, 2015") +
labs(x="", y = "Children given MNP", size=8) +
theme_classic(base_size = 12) + theme(text = element_text(color = "gray20"),
axis.text.x = element_text(NULL),
legend.position = c("bottom"), legend.direction = "horizontal",
legend.justification = 0.05, legend.text = element_text(size = 9, color = "gray10"),
legend.key.height=unit(1,"line"), legend.key.width=unit(1,"line"),
axis.text = element_text(face = "italic"),
axis.title.x = element_text(vjust = -1),
axis.title.y = element_text(vjust = 2),
axis.ticks.y = element_blank(),
axis.line = element_line(color = "gray40", size = 0.5),
axis.line.y = element_blank(),
panel.grid.major = element_line(color = "gray50", size = 0.5),
panel.grid.major.x = element_blank(),
plot.margin = margin(t = 0, r = 0, b = 40, l = 5, unit = "pt"),
plot.title = element_text(face = "bold", color = "black", size = 11, hjust=0.5))

#grid.text("Source: NUTRIDASH 2015", x = .01, y = .03, just = "left", draw = TRUE))
La misma pero interactiva…
library(plotly)
plot_ly(MNP, x = ~Region, y=~childrenMNP, color = ~Region, type = "bar", opacity=0.5)
plot_ly(iycf, x = ~Region, y=~Counselling_mothers, color = ~Region, type = "bar")
Ignoring 115 observationsIgnoring 115 observations
Scatter…
plot_ly(filter(edu_long, area_type=="county"), x = ~coll, y=~HS,
color = ~year, type = "scatter",
text = ~paste("County: ", area_name, 'Year:', year))
Error in filter_(.data, .dots = lazyeval::lazy_dots(...)) :
object 'edu_long' not found
Severe Acute Malnutrition
library(readr)
iycf <- readstata13::read.dta13("data/MAP_2015_General_Info_IYCF.dta")
LS0tDQp0aXRsZTogIkZpbmFsIFByb2plY3QiDQpvdXRwdXQ6IGh0bWxfbm90ZWJvb2sNCi0tLQ0KDQpgYGB7ciwgd2FybmluZz1GQUxTRSwgbWVzc2FnZT1GQUxTRX0NCmxpYnJhcnkobGVhZmxldCkNCmxpYnJhcnkgKHJhc3RlcikNCmxpYnJhcnkocmdkYWwpDQpsaWJyYXJ5KFJDb2xvckJyZXdlcikNCg0KbGlicmFyeShkZXZ0b29scykNCiNkZXZ0b29sczo6aW5zdGFsbF9naXRodWIoInJzdHVkaW8vbGVhZmxldCIpDQoNCnNldHdkKCJDOi9Vc2Vycy9jcmlzaC9PbmVEcml2ZS9Eb2N1bWVudHMvQ29sdW1iaWEgU0lQQS8yMDE3IFNwcmluZy9EYXRhIFZpc3VhbGl6YXRpb24vRmluYWwgUHJvamVjdC8iKQ0KDQpsaWJyYXJ5KHJlYWRyKQ0KaXljZiA8LSByZWFkc3RhdGExMzo6cmVhZC5kdGExMygiZGF0YS9NQVBfMjAxNV9HZW5lcmFsX0luZm9fSVlDRi5kdGEiKQ0Kc2FtIDwtIHJlYWRzdGF0YTEzOjpyZWFkLmR0YTEzKCJkYXRhLzIwMTVfU0FNX1RPVEFMLmR0YSIpDQptbnAgPC0gcmVhZHN0YXRhMTM6OnJlYWQuZHRhMTMoImRhdGEvMjAxNV9NTlAuZHRhIikNCmBgYA0KDQojSW5mYW50IGFuZCBZb3VuZyBDaGlsZCBGZWVkaW5nIChJWUNGKQ0KDQoNCg0KYGBge3J9DQpsaWJyYXJ5KHJhc3RlcikNCndvcmxkX3NwZGYgPC0gc2hhcGVmaWxlKCJkYXRhL3dvcmxkX21hcC9UTV9XT1JMRF9CT1JERVJTX1NJTVBMLTAuMy5zaHAiKQ0KDQojIE1lcmdlIHRoZSBRb0cgZGF0YSB0byBTaGFwZSBGaWxlcw0KbGlicmFyeShkcGx5cikNCklZQ0YyIDwtIHdvcmxkX3NwZGZAZGF0YSAlPiUgbGVmdF9qb2luKGl5Y2YsIGJ5ID0gYyhGSVBTID0gImZpcHMiKSkNCndvcmxkX3NwZGZAZGF0YSA8LSBJWUNGMg0KYGBgDQoNCk1ha2UgbWFwLi4uDQpgYGB7cn0NCmxlYWZsZXQod29ybGRfc3BkZikgJT4lIA0KICAgIHNldFZpZXcobGF0PTEwLCBsbmc9MCAsIHpvb209MikgJT4lDQogICAgYWRkUG9seWdvbnMoDQogICAgc3Ryb2tlID0gRkFMU0UsIGZpbGxPcGFjaXR5ID0gMC41LCBzbW9vdGhGYWN0b3IgPSAwLjUsIA0KICAgIGNvbG9yID0gfmNvbG9yTnVtZXJpYygiUmRZbEduIiwgSGVhbHRoX2ZhY2lsaXRpZXMpKEhlYWx0aF9mYWNpbGl0aWVzKSwNCiAgICBwb3B1cCA9IHBhc3RlKCJDb3VudHJ5OiIsd29ybGRfc3BkZiROQU1FLCI8YnIvPiIsDQogICAgICAgICAgICAgICAgICAiSGVhbHRoIGZhY2lsaXRpZXM6Iixyb3VuZCh3b3JsZF9zcGRmJEhlYWx0aF9mYWNpbGl0aWVzLDEpLCIiLCI8YnIvPiIsDQogICAgICAgICAgICAgICAgICAiQ29tbXVuaXR5IGhlYWx0aCB3b3JrZXJzOiIscm91bmQod29ybGRfc3BkZiRDb21tdW5pdHlfaGVhbHRoX3dvcmtlcnMsMCksIiIpKSAlPiUNCiAgYWRkTGVnZW5kKCJib3R0b21yaWdodCIsIA0KICBwYWwgPSBjb2xvck51bWVyaWMoIlJkWWxHbiIsIHdvcmxkX3NwZGYkSGVhbHRoX2ZhY2lsaXRpZXMpLCB2YWx1ZXMgPSB+SGVhbHRoX2ZhY2lsaXRpZXMsDQogIHRpdGxlID0gIk51bWJlciBvZiBoZWFsdGggY2FyZSBmYWNpbGl0aWVzIiwgb3BhY2l0eSA9IDAuNSkNCg0KYGBgDQoNCg0KIyNNaWNyb251dHJpZW50IHBvd2RlcnMNCg0KIyBDaGlsZHJlbiBnaXZlbiBNaWNyb251dHJpZW50IFBvd2RlciwgcGVyIHJlZ2lvbiwgMjAxNQ0KDQpgYGB7ciwgd2FybmluZz1GQUxTRSwgbWVzc2FnZT1GQUxTRX0NCmNoaWxkcmVuTU5QIDwtIGMoMjAzNTUyLDIyODUyNDIsNTcwNjQ4LDQyNjMyMzQsNDQ2MTk0LDE3ODU4NjcsNjEyMDE2KQ0KUmVnaW9uIDwtIGMoIkNFRS9DSVMiLCAiRWFzdCBBc2lhIGFuZCBQYWNpZmljIiwgIkVhc3Rlcm4gYW5kIFNvdXRoZXJuIEFmcmljYSIsIA0KICAgICAgICAgICAgIkxhdGluIEFtZXJpY2EgYW5kIHRoZSBDYXJpYmJlYW4iLCAiTWlkZGxlIEVhc3QgYW5kIE5vcnRoIEFmcmljYSIsDQogICAgICAgICAgICAiU291dGggQXNpYSIsICJXZXN0IGFuZCBDZW50cmFsIEFmcmljYSIpDQpNTlAgPC0gZGF0YS5mcmFtZShjaGlsZHJlbk1OUCxSZWdpb24pDQpuYW1lcyhNTlApIDwtIGMoImNoaWxkcmVuTU5QIiwiUmVnaW9uIikNCg0KbGlicmFyeShnZ3Bsb3QyKQ0KbGlicmFyeShnZ3RoZW1lcykNCmxpYnJhcnkoc2NhbGVzKQ0KZ2dwbG90KGRhdGE9TU5QLCBhZXMoeD1SZWdpb24sIHk9Y2hpbGRyZW5NTlAsIGZpbGw9UmVnaW9uKSkgKyBnZW9tX2JhcihzdGF0ID0gImlkZW50aXR5IikgKyANCiAgc2NhbGVfZmlsbF9icmV3ZXIocGFsZXR0ZT0iWWxHbkJ1IikgKw0KICBnZ3RpdGxlKCJHcmFwaCAxLiBOdW1iZXIgb2YgY2hpbGRyZW4gZ2l2ZW4gTWljcm9udXRyaWVudCBQb3dkZXIsIHBlciByZWdpb24sIDIwMTUiKSArIA0KICBsYWJzKHg9IiIsIHkgPSAiQ2hpbGRyZW4gZ2l2ZW4gTU5QIiwgc2l6ZT04KSArDQogIHRoZW1lX2NsYXNzaWMoYmFzZV9zaXplID0gMTIpICsgdGhlbWUodGV4dCA9IGVsZW1lbnRfdGV4dChjb2xvciA9ICJncmF5MjAiKSwNCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICBheGlzLnRleHQueCA9IGVsZW1lbnRfdGV4dChOVUxMKSwNCiAgICAgICAgICBsZWdlbmQucG9zaXRpb24gPSBjKCJib3R0b20iKSwgbGVnZW5kLmRpcmVjdGlvbiA9ICJob3Jpem9udGFsIiwNCiAgICAgICAgICBsZWdlbmQuanVzdGlmaWNhdGlvbiA9IDAuMDUsIGxlZ2VuZC50ZXh0ID0gZWxlbWVudF90ZXh0KHNpemUgPSA5LCBjb2xvciA9ICJncmF5MTAiKSwNCiAgICAgICAgICBsZWdlbmQua2V5LmhlaWdodD11bml0KDEsImxpbmUiKSwgbGVnZW5kLmtleS53aWR0aD11bml0KDEsImxpbmUiKSwNCiAgICAgICAgICBheGlzLnRleHQgPSBlbGVtZW50X3RleHQoZmFjZSA9ICJpdGFsaWMiKSwNCiAgICAgICAgICBheGlzLnRpdGxlLnggPSBlbGVtZW50X3RleHQodmp1c3QgPSAtMSksDQogICAgICAgICAgYXhpcy50aXRsZS55ID0gZWxlbWVudF90ZXh0KHZqdXN0ID0gMiksDQogICAgICAgICAgYXhpcy50aWNrcy55ID0gZWxlbWVudF9ibGFuaygpLA0KICAgICAgICAgIGF4aXMubGluZSA9IGVsZW1lbnRfbGluZShjb2xvciA9ICJncmF5NDAiLCBzaXplID0gMC41KSwNCiAgICAgICAgICBheGlzLmxpbmUueSA9IGVsZW1lbnRfYmxhbmsoKSwNCiAgICAgICAgICBwYW5lbC5ncmlkLm1ham9yID0gZWxlbWVudF9saW5lKGNvbG9yID0gImdyYXk1MCIsIHNpemUgPSAwLjUpLA0KICAgICAgICAgIHBhbmVsLmdyaWQubWFqb3IueCA9IGVsZW1lbnRfYmxhbmsoKSwNCiAgICAgICAgICBwbG90Lm1hcmdpbiA9IG1hcmdpbih0ID0gMCwgciA9IDAsIGIgPSA0MCwgbCA9IDUsIHVuaXQgPSAicHQiKSwNCiAgICAgICAgICBwbG90LnRpdGxlID0gZWxlbWVudF90ZXh0KGZhY2UgPSAiYm9sZCIsIGNvbG9yID0gImJsYWNrIiwgc2l6ZSA9IDExLCBoanVzdD0wLjUpKQ0KI2dyaWQudGV4dCgiU291cmNlOiBOVVRSSURBU0ggMjAxNSIsIHggPSAuMDEsIHkgPSAuMDMsIGp1c3QgPSAibGVmdCIsIGRyYXcgPSBUUlVFKSkNCg0KYGBgDQoNCkxhIG1pc21hIHBlcm8gaW50ZXJhY3RpdmEuLi4NCmBgYHtyfQ0KbGlicmFyeShwbG90bHkpDQpwbG90X2x5KE1OUCwgeCA9IH5SZWdpb24sIHk9fmNoaWxkcmVuTU5QLCBjb2xvciA9IH5SZWdpb24sIHR5cGUgPSAiYmFyIiwgb3BhY2l0eT0wLjUpDQpgYGANCg0KYGBge3J9DQpwbG90X2x5KGl5Y2YsIHggPSB+UmVnaW9uLCB5PX5Db3Vuc2VsbGluZ19tb3RoZXJzLCBjb2xvciA9IH5SZWdpb24sIHR5cGUgPSAiYmFyIikNCmBgYA0KDQoNClNjYXR0ZXIuLi4NCmBgYHtyfQ0KcGxvdF9seShmaWx0ZXIoZWR1X2xvbmcsIGFyZWFfdHlwZT09ImNvdW50eSIpLCB4ID0gfmNvbGwsIHk9fkhTLCANCiAgICAgICAgICAgICAgICAgICAgIGNvbG9yID0gfnllYXIsIHR5cGUgPSAic2NhdHRlciIsDQogICAgICAgICAgICAgICAgICAgICB0ZXh0ID0gfnBhc3RlKCJDb3VudHk6ICIsIGFyZWFfbmFtZSwgJ1llYXI6JywgeWVhcikpIA0KYGBgDQoNCg0KI1NldmVyZSBBY3V0ZSBNYWxudXRyaXRpb24NCmBgYHtyfQ0KbGlicmFyeShyZWFkcikNCml5Y2YgPC0gcmVhZHN0YXRhMTM6OnJlYWQuZHRhMTMoImRhdGEvTUFQXzIwMTVfR2VuZXJhbF9JbmZvX0lZQ0YuZHRhIikNCmBgYA0KDQoNCg==